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We report on studies of the structure and dynamics of the (001) surface of single crystal LICU2O2 , 
investigated by He beam scattering at room temperature, and with lattice dynamical models. The 
best fit surface corrugation to measured diffraction patterns shows that the surface termination is 
exclusively a Li +1 Cu +2 -2 plane. Lattice dynamics fits to inelastic He scattering spectra reveal the 
presence of two low-lying surface phonon modes, identified with the motion of Cu +2 , Li +1 surface 
ions normal to the surface. 



I. INTRODUCTION 

LiCu2 02 continues to attract considerable attention 
because of the unique physical properties it exhibits. Ini- 
tially, interest in this system was stimulated by the pres- 
ence of double-chain ladders of Cu +2 0, which presented a 
prototype quasi-one-dimensional (QOD) spin-1/2 quan- 
tum magnetic system with competing magnetic interac- 
tions. Competing magnetic interactions in the double- 
chain ladder were known to give rise to geometric frus- 
tration, which in turn is manifest in an ordered incom- 
mensurate helimagnetic phase at low temperatures. It 
was expected that the presence of S=l/2 spins would 
give rise to strongly competing commensurate quantum 
spin fluctuations that tend to suppress the transition 
temperature^^. More recently, it was discovered that 
this system exhibits ferroelectricity upon the emergence 
of the spiral magnetic order. This renders LiCu202 as 
the second cuprate to join the list of multiferroics^&i. 

A series of successive magnetic phases have been re- 
ported at low temperatures^^: Electron spin resonance 
(ESR) measurements revealed the presence of a dimer- 
ized spin-singlet state at T> 23K, with an energy gap 
of A=72K between the spin-singlet ground state and the 
first spin-triplet excited state. Magnetic neutron scat- 
tering measurements confirmed the existence of a QOD 
spin-ordered helical phase in the temperature range 9K 
<T< 23K 2 ' 8 , and a collinear anti-ferromagnetic (AFM) 
phase was suggested to exist below 9 K 9 ' 10 i n . The pres- 
ence of the QOD magnetic ordering was attributed to 
the presence of impurities in the chains: The observation 
of the classical helical phase in LiCu202 was ascribed to 
substitutional Li + defects that tend to suppress effective 
long-range one-dimensionality^. 

LiCu202 is a quasi-lD insulator. It has a layered 
charge-ordered orthorhombic crystal structure belonging 
to the Pnma space group; the primitive cell has lattice 
constants a=5.73A, 6=2. 86A, c=12.47A respectively^. 
It is a mixed-valent compound with copper ions in the 
Cu +2 and Cu +1 valance states. The magnetic Cu +2 
(S=l/2) ions are located at the center of edge-sharing 
CUO4 plaquettes which form infinite chains along the 
crystallographic 6-axis. Coplanar chains are connected 
by chains of Li ions along the a-axis, and two such planes 




FIG. 1: Crystal structure of LiCu202. The color code is: Li 
(red), O (gray), Cu +1 (blue), Cu+ 2 (green) 



form double-layers parallel to the a6-plane, as shown in 
figure [TJ The QOD spin arrangement is due to these 
double chains of Cu +2 ions that run along the crystallo- 
graphic b axis. The period of each leg of the double spin 
chains is equal to b. The two legs are offset by b/2 rel- 
ative to each other. Along the c-axis, each double-layer 
is separated from its double-layer neighbors by magneti- 
cally inert Cu +1 planes. 

Despite the extensive studies of bulk structural and 
magnetic properties of LiCu202 crystals cited above, no 
investigation of its surface properties has been reported 
in the literature. It is known that LiCu202 crystals easily 
peel along the (001) surface. However, since these crys- 
tals consist of alternating double-layers of Li +1 Cu +2 0^~ 2 
and single layers of Cu +1 , it remains to be determined 
whether the cleavage would result in the coexistence of 
Li +1 Cu +2 02~ 2 and Cu +1 surfaces, or in an exclusive pres- 
ence of one of these layer types. If the latter scenario 
occurs, then it would dictate that the double-layer must 
split in order to provide complete coverage of the two 
newly exposed surfaces. The double-layer splitting means 
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the splitting of the double-chain and the presence of a 
single-chain, which is expected to display different mag- 
netic behavior from that observed in the double-chain. 

In this paper we present the results of extensive studies 
of the (001) surface of LiCu202 , using the experimen- 
tal techniques of clastic and inelastic helium scattering 
aimed at determining its surface structural and dynami- 
cal properties, respectively. Furthermore, we used lattice 
dynamical models with slab geometries to interpret the 
surface phonon dispersion curves derived from the mea- 
sured inelastic scattering spectra. The surface corruga- 
tion topography, derived from a large set of diffraction 
patterns, clearly shows that the surface termination is 
exclusively Li +1 Cu +2 O^T 2 . Moreover, empirical lattice 
dynamics models, with slab geometries based on such 
termination and fit to the inelastic experimental results, 
reveal two low-lying dispersion curves with polarizations 
normal to the surface, one involves the Cu +2 and the 
second involves Li ions. In section II the experimental 
setup and procedures are discussed, and the results and 
conclusion are presented in section III. 



II. EXPERIMENTAL SETUP AND 
PROCEDURE 

Single crystals of LiCu2C>2 with high Li content of 
~ 0.99 ± 0.03 were grown by the floating-zone method. 
The Li content was determined accurately through com- 
bined iodometric titration and thermogravimetric meth- 
ods. This ruled out the possibility of chemical disorder 
between Li and Cu ions. Details of the growth procedures 
and stoichiometry confermation are given in Ref.— 3 . Typ- 
ical crystal samples used were about 3mm x 3mm x 2mm 
in size, with its exposed surface parallel to the afr-plane. 
The crystals were attached to an OFHC copper sample- 
holder by conductive silver epoxy. A cleaving (peeling) 
post was attached to the top sample surface in the similar 
way. The prepared sample holder was mounted on a sam- 
ple manipulator equipped with XYZ motions as well as 
polar and azimuthal rotations. The pressure in the Ultra- 
High Vacuum (UHV) chamber was maintained at 10~ 10 
torr throughout the experiment to ensure the cleanliness 
of the sample surface during measurement performance. 
In situ cleaving under UHV conditions was effected by 
knocking off the cleaving post. Immediately after cleav- 
ing, the quality of the long-range ordering on the sur- 
face was confirmed by the appearance of sharp diffraction 
LEED spots. 

A supersonic mono-energetic collimated helium beam, 
with velocity resolution better than 1.4 %, was generated 
by a nozzle-skimmer assembly and 2mm diameter colli- 
mating slits. The average beam velocity was varied by 
attaching the nozzle reservoir to a closed-cycle helium re- 
frigerator, and controlling the reservoir temperature with 
the aid of a digital temperature controller (Scientific In- 
struments Model 9700) and a diode sensor attached to 
the reservoir. As a result, the beam energy can be var- 



ied in the range 65meV to 21meV by varying the noz- 
zle temperature from 300K to 110K. Polar rotation of 
the sample was used to vary the incident angle 9i with 
respect to the surface normal, while the azimuthal ro- 
tation was employed to align the scattering plane along 
a high-symmetry surface crystallographic direction. The 
scattered He beam was collected by an angle-resolved de- 
tector mounted on a two-axis goniometer, which allows 
the scattered angle Of to be varied independently from 
9r^, and allows in- and out-of the scattering-plane mea- 
surements. The detector— is comprised of an electron 
gun and a multichannel plate (MCP) electron multiplier. 
The electron gun generates a well-collimated, monoener- 
getic electron beam crossing the He beam at right angles. 
The energy of the electron beam is tuned to excite the He 
atoms to their first excited metastable state (2 3 S He*) 
upon impact. Deexcitation of a He* atom at the surface 
of the MCP leads to the ejection of an electron which 
generates an electron cascade that is then collected by 
the anode of the multiplier. By electronically pulsing the 
electron gun, a gate function is created for time-of-flight 
(TOF) measurements in the inelastic HAS mode. The 
details of the detection scheme are given in Ref 14. All 
measurements were performed with the sample surface 
at room temperature. 

By writing the He-atom wave vector as k = (K, k z ), 
where K is the component parallel to the surface, con- 
servation of momentum and energy for in-the-scattering- 
planc geometry can be expressed as 



AK = G + Q = fc/sin(9/ - ki sin 9. 



AE = huj(Q) = E L 



sin 6^ + AK/k, 



sin f 



(1) 
(2) 



where subscripts i and f denote incident and scattered 
beams, respectively, and AK is the momentum trans- 
fer parallel to the surface. G is a surface reciprocal- 
lattice vector, Q is the surface phonon wave vector, 
and Hlu(Q) is the corresponding surface phonon energy. 
Ei — h 2 kf /2M, where M is the mass of a He atom. By 
eliminating k f from the above equations, one obtains the 
so-called scan curve relations which are the locus of all 
the allowed AK and AE as dictated by the conservation 
relations, 



AE = Ei 



sin 6^ + AK/fcj 



smf 



(3) 



The intersections of these scan curves with the phonon 
dispersion curves define the kinematically allowed inelas- 
tic events for a fixed geometric arrangement. Thus, by 
systematically changing E t , 9i, and 9f, the entire disper- 
sion curves can be constructed. 
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The coexistence of both surface terminations would im- 
ply that the intensity of the (1,0) and (0,1) peaks should 
be comparable, but, in fact, the (0,1) peak has exhibited 
appreciably lower intensity than the (1,0) peak. Hence, 
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FIG. 2: Top: Half of a unit cell of LiCu 2 2 . Bottom: Two 
candidates for surface termination: Cu +1 layer with a square 
lattice (bottom right), or a Li +1 Cu +2 OJ 2 layer (bottom left). 
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III. RESULTS AND DISCUSSION 



A. Elastic He Scattering And Surface Structure 

Diffraction patterns were collected from several crystal 
samples at a temperature of 300 K for several scattering 
conditions: 

1. incident He wave numbers, ki, in the rang e 6.42A 
and 11.12A -1 , 

2. incident angles, between 30° and 50°, and 

3. two high-symmetry azimuthal surface orientations, 
a and 6, separated by 90°. 

As was mentioned above, there are two candidates for 
the surface termination, which are shown in Figuref^l the 
Li +1 Cu +2 0^ 2 layer and the Cu +1 layer. The latter has a 
square lattice with lattice constant b, and C*4 V symmetry. 
Thus, the diffraction peaks that correspond to odd mul- 
tiples of 2ir/a should be absent for pure Cu +1 surface ter- 
mination, which contradicts the fact that they do appear 
in all the measured diffraction patterns. Alternatively, 
for Li +1 Cu +2 02~ 2 termination, the well known crystal 
twinningi£i£ would imply that the two azimuthal orien- 
tations would give the same diffraction pattern. This 
means that we obtain a superposition of diffraction pat- 
terns from a(10) and b(01) directions. Moreover, since 
for LiCu2C>2 a ~ 26, the diffraction peaks associated 
with the 6-direction will lie very close to the even-order 
peaks obtained for the a-direction. 

Figure [3] shows a typical diffraction pattern recorded 
for ki = 7.28A" 1 and 9, = 31.4°. Although the twin- 
ning makes the situation a little complicated, luckily, the 
angular resolution in the diffraction pattern allows us to 
clearly separate the peak positions along the two direc- 
tions, as indicated by peaks (2,0) and (0,1) in figured 



FIG. 3: The experimental diffraction pattern with ki = 
7.28A" 1 , 6i = 31.4° at T=300K (red), and the calculated 
diffraction intensities (blue vertical bars) 

The surface periodicity along a and b is derived from 
the positions of their corresponding diffraction peaks us- 
ing the relation 

|ki| (sin^-sin0 G ) = |G| (4) 

where ki is the incident beam wave- vector, G a surface 
reciprocal lattice vector, 9g the corresponding scattering 
angle. Analysis of all the collected diffraction patterns 
yield the surface lattice constant values: a 8 — 5.76 ± 
0.05A and b s = 2.91 ± 0.08A, which compares quite well 
with the bulk values of a=5.73A, b=2.86A and, thus, 
precludes the presence of diffraction satellites and surface 
reconstruction. 

In order to determine the topology of the surface prim- 
itive cell we used the hard corrugated wall model within 
the eikonal approximation to calculate the elastic scat- 
tering intensities 16 . Here, the surface topology is defined 
in terms of a surface corrugation function £(R)— , where 
R denotes a surface position vector and £ the corrugation 
height at that position. In this framework, the scattering 
amplitude of the diffraction peak G is given byi£ 

Ag = ~ f e *[ G ' R +«^( R )] dTL (5) 

b Ju.c. 

where qcz = ki Z — kc z , with kc, z = ki cos#g- The 
integration is carried over the unit cell (u.c.) of area S. 
Since C(R) should have the surface periodicity, we write 

C(R) = E Cgcxp(zG-R) (6) 

G 



We make use of the reflection symmetry perpendicular to 
the 6-direction (j/-axis) of the LiCu202 surface to simplify 
equation ([6]), and write 

((x,y) = ^2 Cn u n 2 sin(ni6ix) cos(n 2 & 2 y) 

+ E C« 3 ,"4 cos(n 3 &ix)cos(n 4 6 2 y) (7) 

n 3 ,n 4 

where b\ = — , b% = =?-, are the surface reciprocal lattice 
vector basis. 

The goal now is to determine the coefficients Cg- The 
following iterative fitting scheme was adopted^ 8 . Ini- 
tially, the magnitudes \Ag\ are determined from the ex- 
perimental diffraction pattern using the geometric rela- 
tion 



P< 



G 



\kcz\ | , ,2 

7i — r a g\ 



(8) 



The experimental intensities are normalized to satisfy the 
unitarity condition 



i 



G 



Next, the eikonal equation 



16,17,18 



A G e jG ' R e 4fcG ^W = 



(9) 



(10) 



G 



where Ac = \Ag\ e^ G , is used to determine C(R). 4>G 
is a diffraction phase angle to be determined. In the first 
iteration we set ^4g = |^-g|, which yields 



Co(R) 



In (-Eg A G e 



iGH\ 



where we replaced fc, by — ki Z in equation (|10p . Further 
iterations involved varying the amplitudes Cg around the 
values obtained from Co(R-)- 

After each iteration the reliability factor R 



R 



cxp 
G 



pcalc\2 
G I 



(ii) 



where N is the number of diffraction peaks, was calcu- 
lated. The iterations were terminated when an accept- 
able small value of R was reached. The corrugation func- 
tion with best-fit parameters, R — 0.003, is 

C(x, y) = -0.06 cos(6ix) - 0.05 cos(2&ix) + 0.02 cos(3&ia;) 
+ 0.03 sin(6ix) - 0.03 cos(b 2 y) + 0.05 cos(26 2 y) 
- 0.07 cos(6 2 y) sin(2&ix) (12) 

It is plotted in figured! and the corresponding calculated 
diffraction peak intensities are shown as blue vertical bars 
in figure [3l The locations of the Li, Cu +2 ,and the two 
O ions have been identified by comparing the positions 
of the corrugation maxima to ionic positions in the unit 
cell; they are indicated in the figure. 
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FIG. 4: Calculated surface corrugation function 
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FIG. 5: Typical time-of-flight spectra. The experimental 
curve (red) is resolved into inelastic peaks (green) and the 
diffuse elastic peak (blue), which is also indicated by arrows. 



B. Inelastic Measurements And Shell-Model 
Calculations 

Inelastic scattering measurements were carried out for 
in-the-scattering-plane geometry. Two typical TOF spec- 
tra showing diffusive elastic as well as inelastic peaks are 
shown in figure [5J The data were collected using beam 
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energies in the range of 25-65 meV. The high-energy 
He beam, Ei = 65 meV, was used to probe the high- 
frequency surface phonon modes. The energy resolution 
for the lowest-energy beam was about 3 meV. Phonon en- 
ergies and momenta were calculated from the TOF peak 
positions, after a Gaussian fit to the peaks, using Eqs. 
(fT|) and @. To characterize the ensuing phonon disper- 
sion points, lattice-dynamical analysis for the bulk and 
surface (slab calculation) was carried out. 



1. Static equilibrium and bulk lattice dynamics 

In lattice-dynamical studies of complex systems such 
as LiCu2C>2 , the construction of a dynamical matrix 
based on pair potentials rather than mere force constants 
is far more advantageous in the sense that the underly- 
ing physics is transparent and many of the potential pa- 
rameters for similar pair interactions can be transferred 
from one compound to another. This is especially useful 
when experimental phonon frequencies are limited to a 
few bulk infrared (IR) and Raman active phonon modes 
at the center of the Brillouin zone (BZ) 14 i 19 , and no 
neutron-scattering data for phonon dispersion is avail- 
able. Bulk lattice dynamics calculations, based on the 
shell-model, was used to obtain an empirical fit to ex- 
isting IR and Raman data. The model incorporates 

TABLE I: Best-fit Born-Meyer and shell parameters. 



Born-Meyer Potential Parameters 



Bond 


a(eV) 


MA- 1 ) 


c(eVA 6 ) 


Cu +1 -0 


18705.6 


2.9 




Cu +2 -0 


10000 


4.4 




Cu +1 -Cu+ 2 


3490.95 


3.0 




Cu +1 -Cu +1 


18853 


3.0 




Cu +2 -Cu+ 2 


6598.44 


3.0 




o-o 


2146.30 


3.4 


50 


Li-O a 


840 


3.0 





shell model parameters 



Ion 


Z(e) 


Y(e) 


K(eV/A 2 ) 


Li 


+1 


1.5 


5 


Cu +1 


+1 


1.7 


20 


Cu+ 2 


+2 


2.9 


30 


o b 


-2 


-3.1 


40 



refiV ref.il 



two-body central potentials, namely the Coulomb poten- 
tial V$ (r) = ZiZje 2 /r for the long-range interactions 
and either Born-Mayer Vij(r) = a,ije~ bijr or Bucking- 
ham Vij(r) = a,ij exp(-bijr) — Cij/r 6 potentials for the 
short-range interactions. The ionic shells are character- 
ized by the ionic charge Z, the shell charge Y, and the 
intra-ion shell-core force constant K. Some of the initial 
short-range pair potentials were obtained from the litera- 
ture: (Cu +2 -0) 14 , (O-O) 14 and (Li-O) 19 . The remaining 
pair potential parameters were set to satisfy the static 
equilibrium conditions 2 - which state that the forces on 
the particles in their equilibrium positions should van- 
ish. This treatment also ensures the consistency of the 
static and dynamical properties of the crystal. It should 

TABLE II: Comparison of experimental and calculated IR 
and Raman frequencies 



Raman active 



Mode 


Polarization 


Frequency (cm ) 






Experiment- 


Calculated 






573 


593.30 






497 








460 


476.76 


A g 


aa 


367 


375 






297 


275.9 






178 


191.33 






167 








122 


108.86 



IR active 



Mode 


Polarization 


Frequenc 
Experiment^ 


y (cm *) 
Calculated 


B2u 


b 


240 
288 
312 
416 


268.79 
449.88 


B3 U 


a 


240 
288 
320 
392 
440 


222.88 
269.621 

375.54 
463.08/474.46 



be noted, however, that satisfying static equilibrium con- 
ditions does not guarantee dynamical stability, namely 
ensuring the reality of the phonon frequencies through- 
out the BZ. It should also be noted that the shell-model 
parameters do not appear in the static equilibrium equa- 
tions ions are treated as rigid bodies. However, these pa- 
rameters are introduced into the dynamical matrix and 
are determined through the process of fitting the 17 ex- 



6 



perimental IR and Raman modes reported in the litera- 
ture. The short-range potential and shell-model parame- 
ters that produce the best fit for these modes are listed in 
Table I. A comparison of the experimental and calculated 
values of these modes is given in Table II. 

The complete set of calculated bulk phonon dispersion 
curves (60 branches) is shown in figure along the high- 
symmetry directions T-X , T-Y and T-Z. 

According to the point group mmm (D 2 h) of the space 
group Prima, there are 30 Raman active and 27 IR active 
modes at the T-point, which are classified in terms of the 
corresponding irreducible representations as 

10A g © 5B lg © 10B 2g © 5B 3g Raman 
9B lu © 4B 2u © 9B 3u IR 

The remaining modes at the T-point include five A u 
modes and 3 zero-frequency acoustic modes. The identi- 
fication of the symmetries of the calculated modes at the 
T-point was established with the aid of the symmetry 
projection operators 2 ^ of the point group mmm. Away 
from the T-point, the point group of the wave vector is 
isomorphic with C 2v for all 3 high-symmetry directions. 
Symmetry compatibility relations give 



c 2v 


D 2h 




A ((100)) 


£((010)) 


A ((001)) 


Ai 


Ag 


Ag 

B 2u 


Ag 

Bi„ 


A 2 


A M 
B 3g 


A u 

B 2 g 


A M 

Big 


Bi 


Big 

B 2u 


Big 

^3u 


B 3 g 

B 2ll 


B 2 


B 2 g 

Bi„ 


B 3 g 

Biu 


B 2 g 

B3 M 



In all, we have 20Ai © 10A 2 © 10Bi © 20B 2 along the 
(100)- and the (OOl)-directions, and 15Ai©15A 2 ©15Bi© 
15B 2 along the (OlO)-direction. 

The highest two bands (~ 800-1200cm" 1 ) are 
longitudinal optic phonons involving the motion 
of the eight 0~ 2 ions in the primitive cell; in 
the upper band the four 01 ions at positions 
{0.137,1/4,0.405}, {0.363,3/4,0.905}, {0.863,3/4,0.595}, 
{0.637,1/4,0.095} move in-phase with each other, 
but anti-phase with the four 02 ions at positions 
{0.115, 1/4, 0.105}, {0.385, 3/4, 0.605}, {0.885, 3/4, 0.895}, 
{0.615, 1/4,0.395}, while in the lower band, the 01 and 
02 ions still have anti-phase motion, but also the 
ions at y — 1/4 have anti-phase motion with those at 
y = 3/4. The lowest three bands (acoustic phonons) 



mix with higher bands at higher q vectors in the a- and 
6-directions, but remain distinct in the c-direction. 



2. Lattice Dynamics of Slab-geometry and 
Surface Phonons 

A lattice dynamics shell-model with slab-geometry con- 
sisting of 36 layers, and terminated with Li +1 Cu +2 0^~ 2 
surfaces, was used to calculate the phonon dispersion 
curves of that surface. Since translation symmetry is 
broken normal to the slab surfaces, primitive and non- 
primitive translations along this direction are not al- 
lowed. Consequently, the slab symmetry reduces to 2ma, 
while the surface symmetry is plm. 

Surface equilibrium analysis showed that the Cu +2 and 
Li +1 surface ions have to be displaced outward, along the 
surface normal, from its bulk position, in order to satisfy 
equilibrium conditions in that direction. The pair po- 
tential and shell-model parameters of Table I were used 
in the initial slab calculations. However, the surface pa- 
rameters had to be modified in order to achieve global 
stability in the entire surface BZ. 

TABLE III: Character of the surface phonon dispersion 
curves. (Pol.:polarization, AP: anti-phase, IP: in-phase) 



Branch 


(10) 


(01) 




Type 




Ions 


Pol. 


Ions 


Pol. 


Si 


Li +1 


2 


Cu+ 2 , Li +1 
1 : 5 


z 


AP ((01)) 


s 2 


Cu+ 2 


Z 


Cu+ 2 , Li +1 
5 : 1 


z 


IP ((01)) 


s 3 


Cu+ 2 


y 


Cu+ 2 


X 




S* 


Cu+ 2 


X 


Cu+ 2 


y 




Ss 


o- 2 


z 


o- 2 


z 




s 6 


Li +1 


y 


Li +1 


X 




s 7 


Li +1 


X 


Li +1 


y 




s 8 


o- 2 


y 


o- 2 


X 


IP 


S 9 


o- 2 


y 


o- 2 


y 


AP 


Sio 


o- 2 


X 


o- 2 


X 


AP 



The resulting surface phonon dispersion curves are 
shown as blue dots superimposed on a gray background 
in figure [2 they are labeled Si, i = 1, . . . , 10. The gray 
background areas and lines resulting from the slab cal- 
culations correspond to contributions from the slab bulk, 
namely projections onto the surface BZ. Some of the dis- 
persed gray lines would turn into solid bands as the thick- 
ness of the slab approaches infinity. The general charac- 
ter of the surface modes Si — Sio is given in Table III. Si 
and S 2 involve the motions of the Cu +2 and Li +1 normal 
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to the surface. While S\ is quite flat for wave vectors 
q ~ 0.15 — 0.5(27r/a) along the a-direction, it exhibits 
anomalous softening in the proximity of the surface BZ 
boundary along the 6-direction. This softening could be 
the result of the tenuous equilibrium of these ions normal 
to the surface. 

The scattering of He beams with energies in the range 
25-65 meV allows measurement of surface phonons with 
frequencies bellow 400 cm -1 (50 meV). Thus we will 
focus on identifying the measured inelastic scattering 
events with calculated surface phonon dispersion curves 
in that range. The process of identification is further 
complicated by the presence of twining in the sample 
crystals. A succesful procedure must be capable of sort- 
ing out modes propagating in the a-dircction from those 
propagating in the 6-direction. Figure 8(a) illustrates this 
complexity Here, all measured inelastic events are plot- 
ted in the left and right panels, after reducing their mo- 
menta to the first surface BZ along the a-direction (left 
panel) and along the 6-direction (right panel). Notice 
that each measured inelastic event has a different value 
of its reduced momentum for the different directions. The 
sorting criterion adopted, was to assign a given event to 
the direction where it is nearest in energy to a disper- 



sion curve. The final result is plotted in figure 8(b). The 
agreement is surprisingly quite good, especially with the 
two low-lying dispersion curves. In addition, this agree- 
ment supports the surface termination analyzed in this 
paper. 

In summary, with the aid of He scattering techniques 
at room temperature, we find that the (001) surface of 
LiCu2C>2 is exclusively terminated by Li +1 Cu +2 0^ 2 , and 
that no surface reconstruction occurs. Empirical fit- 
ting of lattice dynamics shell-model based calculations 
to measured inelastic He scattering spectra, supports the 
proposed termination and reveals that the lowest surface 
phonon dispersion branches involve the motion of Cu+ 2 
and Li +1 ions normal to the surface. 
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FIG. 6: Bulk phonon dispersion curves along a (£ = q/G a ), b (£ = q/Gb) and c (£ = q/G c )- 
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FIG. 8: Top panels: All measured inelastic events, reduced 
to the proper surface BZ, superposed on calculated dispersion 
curves along (10), (left panel), and (01) (right panel). Lower 
panels: Fits obtained by the procedure described in the main 
text for (10), (left panel), and (01) (right panel). (Error bars 
are also indicated in the lower panels). 



